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Abstract. In recent years quantum statistical mechanics have benefited of 
cultural interchanges with quantum information science. There is a bulk of 
evidence that quantifying the entanglement allows a fine analysis of many relevant 
properties of many-body quantum systems. Here we review the relation between 
entanglement and the various type of magnetic order occurring in interacting spin 
systems. 
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1. Introduction 

In one of the most influential papers for the foundation of quantum 
physics (Schrodinger 1935), entanglement was recognized as not "one but rather the 
characteristic trait of quantum mechanics" . Although such kind of non-local corre- 
lations were thoroughly explored for the analysis of the conceptual fundaments of 
quantum mechanics (Bell 1987, Peres 1993), the interest in understanding the proper- 
ties of entangled states has received an impressive boost with the advent of quantum 
information. In fact it was understood that non-local correlations are responsible for 
the enhanced efficiency of quantum protocols (Nielsen & Chuang 2000) . In some other 
cases, quantum teleportation just to mention an important example, quantum corre- 
lations are a necessary ingredient. Entanglement can thus be considered a resource 
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for quantum information processing. The need for denning the amount of resources 
required to implement a given protocol has lead to a fertile line of research in quantum 
information aiming at the quantification of entanglement (Horodecki et al. 2007). To 
this end, necessary criteria for any entanglement measure to be fulfilled have been 
elaborated and lead to the notion of an entanglement monotone (Vidal 2000). There 
is a fairly clear scenario in the case of bipartite systems. The multipartite counterpart 
is much more intricate. 

The large body of knowledge developed in the characterization of entanglement 
has found important applications in areas that, less than a decade ago, were quite 
distant from quantum information. This is the case of quantum statistical mechanics. 
Traditionally the characterization of many-body systems has been carried on through 
the study of different physical quantities (as the magnetization in magnetic systems) 
and their correlations. Very little attention was paid to the structure of their quantum 
state and in particular to its amount of entanglement. The research at the border 
between these areas is rapidly evolving and has lead to numerous interesting results. 
The promise of this new interdisciplinary field is twofold: firstly, to improve our 
understanding of strongly correlated systems, beyond the current state of the art; 
and secondly, to control and manipulate quantum correlations to our convenience, 
ultimately with the hope of making an impact on the computation schemes to solve 
'hard problems' in computer science. 

Methods developed in quantum information have proven to be extremely useful 
in the analysis of the state of many-body systems. At the same time experience 
built up over the years in condensed matter is helping in finding new protocols 
for quantum computation and communication. The cross-fertilization of quantum 
information with statistical mechanics should not, however, come as a surprise. 
Quantum computers are themselves many-body systems, the main difference from 
traditional solid-state systems is that a quantum computer can be controlled and 
operate under non-equilibrium conditions. It is therefore natural to profit of the 
best instruments developed so far in the two disciplines for the understanding of 
quantum complex systems, being quantum computers or condensed matter systems. 
In this review we deal with a particular aspect of this research area: We will focus on 
interacting spin systems on a lattice and describe the relations between entanglement 
and magnetic ordering. A more detailed discussion on these issues can be found 
in (Amico et al. 2008, Vedral 2008). 

A complete characterization of entanglement in a many-body system is hopeless. 
The number of possible ways in which the system can be partitioned explodes 
on increasing the number of elementary constituents (namely spins). Nevertheless 
judicious choices have made possible to highlight interesting properties of many-body 
correlations. Numerous contributions in this volume consider a bipartition in which 
the system is divided into two distinct regions. If the total system is in a pure state, as 
when the system is in its ground state, then a measure of the entanglement between the 
two regions is given by the von Neumann entropy associated to the reduced density 
matrix (of one of the two regions). This is not however the only possible choice, 
our paper will briefly review other way to quantify entanglement whose properties 
contribute to our understanding of many-body systems. 

In the case of bipartite entanglement one can consider the quantum correlation 
between two given spins after having traced out the rest of the system. In this 
case the entanglement between the two selected sites can be quantified by the 
concurrence (Wootters 2001). The study of two-site entanglement, as we will briefly 



CONTENTS 



3 



describe in the following of this paper, allows to detect quite well the presence of 
quantum phase transition in the phase diagram. A different approach to entanglement 
in many-body systems arises from the quest to swap or transmute different types of 
multipartite entanglement into pairwise entanglement between two parties by means 
of generalized measures on the rest of the system. In a system of interacting spins on 
a lattice this means to maximize the entanglement between two spins by performing 
measurements on all the others. With this aim, the concept of localizable entanglement 
has been introduced in (Verstraete, Martin-Delgado & Cirac 2004, Popp et al. 2005). 

The structure of a many-body system's quantum state is by far much richer 
than that captured by bipartite entanglement. In the multipartite case the grounds 
for quantitative predictions are less firm because of the exceptional difficulty of 
the problem. Nevertheless there are a number of very interesting results already 
available. Among all we mention here the bounds which have been derived on 
the ground state energy (Giihne et al. 2005) which allow to discriminate among 
different n-particle quantum correlations. The idea of deriving bounds for macroscopic 
quantities (as for example the ground state energy) is a quite powerful methods 
to characterize entanglement in many-body system and is related to the concept of 
entanglement witness. As an interesting connection between statistical mechanics and 
quantum information it turns out that in many cases thermodynamic quantities as 
the magnetization or the susceptibility (in the case of interacting spin model) behave 
as entanglement witness (Toth 2005, Wiesniak et al. 2005) thus providing a way to 
detect entanglement experimentally. Finally we would like to mention that the natural 
dynamics of condensed matter systems may be important to detect bound entangled 
states (Patane et al. 2007, Ferraro, Cavalcanti, Garcia-Saez & Acin 2008) which are, 
in several cases, difficult to be realized artificially. 

Topic of this short review is to complement the different papers in this volume 
by describing various measure of entanglement, other then the block entropy, which 
were used to characterize the equilibrium and dynamical properties of entanglement 
in spin systems (see also the article by J. I. Latorre and A. Riera in the same volume). 
Interacting spin models (Auerbach 1998, Schollwock et al. 2004) provide a paradigm 
to describe a wide range of many-body systems. They account for the effective 
interactions in a variety of very different physical contexts ranging from high energy 
to nuclear physics. In condensed matter beside describing the properties of magnetic 
compounds, they capture several aspects of high-temperature superconductors, 
quantum Hall systems, heavy fermions, just to mention few important examples. Of 
course interacting spins are central to quantum information processing (Nielsen & 
Chuang 2000). 

The paper is organized as follows. In the next Section we list the spin models 
that will be considered thereafter. Then we give a brief overview of the entanglement 
measures that are currently used to analyze many-body systems. In the Section 4 we 
discuss the main outcomes of the uses of this entanglement measure applied to the 
spin models. The conclusions and a possible outlook will be presented in Section 5. 
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2. Model systems 

A large class of relevant models of interacting spin on a e?-dimensional lattice can be 
described by the Hamiltonian 

h=1 J2 [ J x j)s i s j + 4 ij)s i s j + 4 ij)s i s j] - hZ Yl s " ■ w 

In Eq.(l) Sf (a = x,y,z) are spin-1/2 operators defined on the i-th site of a d- 
dimensional lattice. The ground state of Eq.(l) is, in general, a highly entangled 
state. Nonetheless it exists a point where the ground state is indeed a classical state, 
factorizcd in the direct space (Kurmann et al. 1982, Roscilde et al. 20056, Giampaolo 
et al. 2009). Such phenomenon occurs at a precise value of the field hf that 
was obtained in any dimensional bipartite lattice and for finite range exchange 
interaction (Giampaolo et al. 2009) and even in presence of frustration (Giampaolo 
et al. 2009). Despite the similarities with the saturation phenomenon occuring 
in ferromagnets in external magnetic fields, it was proved that the factorization 
of the ground state is due to a fine tuning of the control parameters, within the 
magnetically ordered phase. Interestingly enough, the results obtained so far indicate 
that the factorization point is a precursor of the quantum phase transition. We 
also mention that analysis of the ground state for finite size systems demonstrated 
that the factorization can be viewed as transition between ground states of different 
parity(Rossignoli, Canosa & Matera 2008, Giorgi 2009). 

For nearest neighbors, the previous Hamiltonian defines the XYZ anisotropic 
Heisenberg model. In this case the exchange couplings are commonly parameterized 
as J x = J(l + 7), J y = J(l — 7), and J z = 2JA. A positive (negative) exchange 
coupling J favors antiferromagnetic (ferromagnetic). In one dimension the model 
defined by the Hamiltonian 1 is exactly solvable in several important cases. This is 
particularly interesting in the analysis of entanglement because quantum effects are 
particularly pronounced in low dimensions. 

Whenever A — the (quantum anisotropic XY) Hamiltonian can be diagonalized 
by first applying the Jordan- Wigner transformation and then performing a Bogoliubov 
transformation (Lieb et al. 1961, Pfeuty 1970, Barouch & McCoy 1971). The quantum 
Ising model corresponds to 7 = 1 while the (isotropic) AA-model is recovered for 
7 = 0. In the isotropic case the model possesses an additional symmetry resulting 
in the conservation of the total magnetization along the z-axis. The properties of 
the Hamiltonian are governed by the dimensionless coupling constant A = J/2h. The 
phase diagram is sketched in Fig. 1 In the interval < 7 < 1 the system undergoes 
a second order quantum phase transition at the critical value A c = 1. The order 
parameter is the magnetization in x-direction, (S x ), different from zero for A > 1. In 
the phase with broken symmetry the ground state has a two-fold degeneracy reflecting 
a global phase flip symmetry of the system. The magnetization along the z-direction, 
(S z ), is different from zero for any value of A, but is singular behavior in its first 
derivative at the transition point. In the whole interval < 7 < 1 the transition 
belongs to the Ising universality class. For 7 = the quantum phase transition is of the 
Berezinskii-Kosterlitz-Thouless type. In several cases the evaluation of entanglement 
requires the determination of the average magnetization Mf(t) = {tjj\Sf {t)\^p) and 
of the equal-time correlation functions g"^(t) = ( , ip\Si'(t)S^ l (t)\il)). These correlators 
have been calculated for this class of models in the case of thermal equilibrium (Lieb 
et al. 1961, Pfeuty 1970, Barouch & McCoy 1971). These can be recast in the form of 
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Figure 1. The zero temperature phase diagram of the one dimensional 
anisotropic XY model in transverse field. Along the quantum critical line the 
model identifies the Ising universality class with indices z = v = 1; in the hatched 
area (with 7 > 0) the system display long range order in the x—y spin components. 
The critical XY regime coicide with that one of the XXZ model for A = 0. The 
factorization of the ground state occurs along the circle A -1 = y/l — 7 2 . 




Figure 2. Zero temperature phase diagram of the spin 1/2 XXZ model in 
one dimension. The XY phase is characterized by power law decay of the xy 
correlations. The Neel and the XY phases are separated by a line of second order 
phase transitions; at A = 1 the transition is of Berezinskii-Kosterlitz-Thoulcss. 
The ferromagnetic and XY phases are separated by first order phase transitions 
due to simple level crossings; the onset to the ferromagnetic phase occurs through 
the saturation phenomenon. 



Toeplitz determinants in equilibrium case and can be expressed as a sum of PfafHans 
in certain non-equilibrium situations (Amico & Osterloh 2004). 

In the case in which 7 = and for any value of A the model is referred to 
as the XXZ model. The two isotropic points A = 1 and A = — 1 describe the 
antifcrromagnetic and ferromagnetic chains respectively (see the phase diagram in 
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Fig. 2). The one-dimensional XXZ Heisenberg model can be solved exactly by Bethe 
Ansatz technique (see e.g. (Takahashi 1999)) and the correlation functions can be 
expressed in terms of certain determinants (see (Bogoliubov et al. 1993) for a review). 
Correlation functions, especially for intermediate distances, are in general difficult 
to evaluate, although important steps in this direction have been made (Kitaninc 
et al. 1999, Gohmann & Korepin 2000, Boos et al. 2008). The zero temperature 
phase diagram of the XXZ model in zero magnetic field shows a gapless phase in 
the interval — 1 < A < 1 with power law decaying correlation functions (Mikeska 
& Pesch 1977, Tonegawa 1981). Outside this interval the excitations are gapped. 
The two phases are separated by a Berezinskii-Kosterlitz-Thouless phase transition at 
A = 1 while at A = — 1 the transition is of the first order. In the presence of the 
external magnetic field a finite energy gap appears in the spectrum. The universality 
class of the transition is not affected, as a result of the conservation of the total spin-z 
component (Takahashi 1999). 

Another interesting case of the Hamiltonian in Eq.(l) is when each spin interacts 
with all the other spins in the system with the same coupling strength 

ij i 

For site-independent magnetic field hf = h a \fi, a — x,y,z, this model is known 
as the Lipkin-Meshkov-Glick (LMG) model (Lipkin et al. 1965, Meshkov, Lipkin 
& Glick 1965, Meshkov, Glick & Lipkin 1965). In this case the dynamics of the 
system can be described in terms of a collective spin S a — J2j At J 2 1 = 4/i z 
the Hamiltonian manifests a supersimmetry (Unanyan & Fleischhauer 2003). The 
phase diagram depends on the parameter A = J/2h z . For a ferromagnetic coupling 
(J > 0) and h x = h v = the system undergoes a second order quantum phase 
transition at A c — 1, characterized by mean field critical indices (Botet et al. 1982). 
For h v — 0, h z < 1 and 7 = the model exhibits a first order transition at 
h x = (Vidal et al. 2006) while for an antiferromagnetic coupling and h v ~ a first 
order phase transition at h z — occurs for any 7's. The model Hamiltonian defined 
above embraces an important class of interacting fermion systems with pairing force 
interaction (like the BCS model). Both the LMG and the BCS type models can be 
solved exactly by Bethe Ansatz (Richardson 1963, Richardson & Sherman 1964). 

We end this very brief overview with spin-1 systems which where originally 
considered to study the quantum dynamics of magnetic solitons in antifcrromagnets 
with single ion anisotropy (Mikeska 1995). In one dimension, half-integer and integer 
spin chains have very different properties (Haldanc 1983a, Haldanc 19836). We will 
see that the typical ground state of such models displays characteristic features in its 
entanglement content. Long range order that is established in the ground state of 
systems with half- integer spin (Lieb et al. 1961), may be washed out for integer spins. 
In this latter case, the system has a gap in the excitation spectrum. A paradigm 
model of interacting spin-1 systems is 

JV 

ff = ]Ts j .S m +/3(S J .S m ) 2 (2) 

i=0 

The lack of long range order arises because of the presence of zero as an eigenvalue of 
Sf; the corresponding eigenstates represent a spin configuration that can move freely 
in the chain, ultimately disordering the ground state of the system, experiencing a 
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gap with the lowest energy excitation (Mikeska 1995, Gomez-Santos 1991). The so 
called string order parameter was proposed to capture the resulting 'floating' Neel 
order, made of alternating spins | f), | J.) with strings of |0)'s in between (den Nijs & 
Rommelse 1989) 

i+R-l 

0% ring =limJS?( J] e™ s %)S? +R ) . (3) 

fe=i+l 

The ground state of physical systems described by Hamiltonians of the form of Eq.(2) 
has been studied in great details (Schollwock et al. 2004). Different phase transitions 
have been found between antiferromagnetic phases, Haldane phases, and a phase 
characterized by a large density of vanishing weights (Sf — 0) along the chain. Some 
features of the phenomenology leading to the destruction of the antiferromagnetic 
order can be put on a firm ground for [3 = 1/3 (AKLT model), where the ground state 
of the Hamiltonian in Eq.(2) is known exactly (Affleck et al. 1988). In this case it was 
proved that the ground state is constituted by a sea of nearest neighbor valence bond 
states, separated from the first excitation by a finite gap with exponentially decaying 
correlation functions. 

3. Entanglement measures 

The study of quantum correlations in many-body systems depends heavily on 
the impressive progress that has been achieved in the theory of entanglement 
quantification. The new ingredient that makes the many-body case very appealing 
is the rich variety of ways in which the system can be partitioned into subsystems. 
Comprehensive overviews of entanglement measures can be found in (Bcngtsson 
& Zyczkowski 2006, Brufi 2002, Eiscrt 2006, Horodccki et al. 2007, Plenio & 
Vedral 1998, Plenio & Virmani 2007, Vedral 2002, Wootters 2001), it is however 
convenient to recall some of the entanglement measures that are routinely used 
to characterize many-body systems. Important requirements for an entanglement 
measure are that it should be invariant under local unitary operations; it should be 
continuous and, furthermore, additive when several identical copies are considered. 

Most of the work done on entanglement in many-body systems deals with the 
bipartite case. A pure bipartite state is not entangled if and only if it can be 
written as a tensor product of pure states of the parts. It can be demonstrated that 
reduced density matrices can be decomposed by exploiting the Schimdt decomposition: 
Pb/a — a i | ^pB/A.i) {^b/a^ | ■ Since only product states lead to pure reduced 
density matrices, a measure for their mixedness points a way towards quantifying 
entanglement. One can thus use a suitable function of the a, given by the Schmidt 
decomposition to quantify the entanglement. Remarkably enough the von Neumann 
entropy S(pb/a) = J2i a i log(ctf ) , can quantify the entanglement encoded in Pb/a- 
We point out that an infinite class of entanglement measures can be constructed for 
pure states. Infact by tracing out one of two qubits in the state, the corresponding 
reduced density matrix pa contains only a single independent parameter: its 
eigenvalue < 1/2. This implies that each monotonic function [0,1/2] i— > [0,1] of 
this eigenvalue can be used as an entanglement measure. A relevant example is the 
(one-) tangle (Coffman et al. 2000) ti[pa] = 4:dctpA- By expressing pa in terms 
of spin expectation values, it follows that t^/j^] = \ — ((S x ) 2 + (S v ) 2 + (S z ) 2 ) 
where (S a ) = tr A {p A S a ) and S a = \a a , a a {a = x,y,z} being the Pauli 
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matrices, For a pure state of two qubits it can be shown that n is equivalent to 
the concurrence C (Hill & Wootters 1997, Wootters 1998) for pure states of two 
qubits. The von Neumann entropy can be expressed as a function of the (one-) tangle 
S[pa] = h[(l + y/l — ti{pa])/2] where h(x) —: — xlog 2 x — (1 — x) log 2 (l — x) is the 
binary entropy. 

Subsystems of a many-body (pure) state will generally be in a mixed state. 
In this case different ways of quantifying entanglement can be realized. Three 
important representatives are the entanglement cost Eq, the distillable entanglement 
Ed (both defined in Rcf. (Bennett, Bernstein, Popescu & Schumacher 1996)) and 
the entanglement of formation Ep (Bennett, DiVincenzo, Smolin & Wootters 1996). 
In the following we concentrate on the entanglement of formation. The conceptual 
difficulty behind its calculation lies in the infinite number of possible decompositions 
of a density matrix. Therefore, even knowing how to quantify bipartite entanglement 
in pure states, we cannot simply apply this knowledge to mixed states in terms of an 
average over the mixtures of pure state entanglement. It turns out that the correct 
procedure is to take the minimum over all possible decompositions. This conclusion 
can be drawn from the requirement that entanglement must not increase on average 
by means of local operations including classical communication. The entanglement of 
formation of a state p is therefore defined as 

E F (p) := mm^2 Pj S(p Atj ) , (4) 

3 

where the minimum is taken over all realizations of the state pab = Pj|^j)(V ; j|> 
and S(paj) is the von Neumann entropy of the reduced density matrix paj ■= 
tr b | i>j) (4>j |- For systems of two qubits, an analytic expression for Ep does exist 
and it is given by 

E FM = -Z^^^^^ <») 

<T = ± 

where C(p) is the the so called concurrence (Wootters 1998, Wootters 2001) defined 

as 

C = max{Ai - A 2 - A 3 - A 4 , 0} . (6) 

where A? > . . . > A4 are the eigenvalues of R = ^fpp^fp = ^/p(a v <Sia y )p*(a v ®a y )^/p (* 
indicates the complex conjugation). As the entanglement of formation is a monotonous 
function of the concurrence, also C itself or its square T2 - called also the 2-tangle - 
can be used as entanglement measures. The concurrence C and the tangle n both 
range from (no entanglement) to 1. By virtue of (6), the concurrence in a spin-1/2 
chain can be computed in terms of the two-point spin correlation functions. As an 
example (that is relevant for the present article) we consider a case where the model 
has a parity symmetry, it is translational invariant and the Hamiltonian is real. In 
this case the concurrence reads 

dj : =2max{0,C4-,C4 / } . (7) 

where C& = \gf? + g™\ - ^(1/4 + , 9 f/) 2 - M* and C// = \gf? - g™\ + gff - 1/4, 
with gff = (S?S?) and M z = (S z ). A state with dominant fidelity of parallel and 
anti-parallel Bell states is characterized by dominant C 1 and C 11 , respectively. This 
was shown in (Fubini et al. 2006), where the concurrence was expressed in terms of the 
fully entangled fraction as defined in (Bennett, DiVincenzo, Smolin & Wootters 1996). 
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The importance of the tangle and the concurrence is due to the monogamy of 
entanglement which was expressed in terms of an inequality in (Coffman et al. 2000) 
for the case of three qubits. This inequality has been proved to hold also for n- 
qubits system (Osborne & Verstraete 2006). In the case of many-qubits it reads 
J2j^i^ij — r M ■ The so called residual tangle n,, — ^jj^Cfj, is a measure for 
multipartite entanglement not stored in pairs of qubits only. 

Another measure of entanglement we mention is the relative entropy of 
entanglement (Vedral et al. 1997). It can be applied to any number of qubits in 
principle (or any dimension of the local Hilbcrt space). It is formally defined as 
E(<t) := min pe x) S(<r\\p), where S^ctH/?) = tr er[mcr — \np] is the quantum relative 
entropy. This relative entropy of entanglement quantifies the entanglement in a by 
its distance from the set V of separable states. The main difficulty in computing 
this measure is to find the disentangled state closest to p. This is in general a non 
trivial task, even for two qubits. In the presence of certain symmetries - which is 
the case for e.g. eigenstates of certain models - an analytical access is possible. In 
these cases, the relative entropy of entanglement becomes a very useful tool. The 
relative entropy reduces to the entanglement entropy in the case of pure bi-partitc 
states; this also means that its, so called, convex roof extension (Uhlmann 1998) 
coincides with the entanglement of formation, and is readily deduced from the 
concurrence (Wootters 1998). 

It is important to realize that not just the quantification of many-party 
entanglement is a difficult task; it is an open problem to tell in general, whether a state 
of n parties is separable or not. It is therefore of great value to have a tool that is able to 
merely certify if a certain state is entangled. An entanglement witness W is a operator 
that is able to detect entanglement in a state. The basic idea is that the expectation 
value of the witness W for the state p under consideration exceeds certain bounds only 
when p is entangled. An expectation value of W within this bound however does not 
guarantee that the state is separable. Nonetheless, this is a very appealing method 
also from an experimental point of view, since it is sometimes possible to relate the 
presence of the entanglement to the measurement of few observables. Simple geometric 
ideas help to explain the witness operator W at work. Let T be the set of all density 
matrices and let £ and T> be the subsets of entangled and separable states, respectively. 
The convexity of I? is a key property for witnessing entanglement The entanglement 
witness is then an operator defining a hyper-plane which separates a given entangled 
state from the set of separable states. The main scope of this geometric approach is 
then to optimize the witness operator (Lewenstein et al. 2000) or to replace the hyper- 
plane by a curved manifold, tangent to the set of separable states (Giihne 2004). We 
have the freedom to choose W such that tr (puW) < for all disentangled states 
Pd G T>. Then, tr (pW) > implies that p is entangled. Entanglement witnesses are 
a special case of a more general concept, namely that of positive maps. These are 
injective superoperators on the subset of positive operators. When we now think of 
super-operators acting non-trivially only on part of the system (on operators that act 
non trivially only on a sub-Hilbert space), then we may ask the question whether a 
positive map on the subspace is also positive when acting on the whole space. Maps 
that remain positive also on the extended space are called completely positive maps. 
Positive but not completely positive maps are important for entanglement theory. 
Indeed it can be shown (Horodecki et al. 1996) that state pab is entangled if and only 
if a positive map A exists (not completely positive) such that (1U <g> Nb)pab < 0. For 
a two dimensional local Hilbert space the situation simplifies and in a system of two 
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qubits the lack of complete positivity in a positive map is due to a partial transposition. 
This partial transposition clearly leads to a positive operator if the state is a tensor 
product of the parts. In fact, also the opposite is true: a state of two qubits p AB is 
separable if and only if p\\ > that is, its partial transposition is positive. This is very 
simple to test and it is known as the Peres-Horodecki criterion (Peres 1996, Horodecki 
et al. 1996). The properties of entangled states under partial transposition lead to a 
measure of entanglement known as the negativity. The negativity Nab of a bipartite 
state is defined as the absolute value of the sum of the negative eigenvalues of p T A A B . 
The logarithmic negativity is then defined as 

E N = log 2 2(2iV Ai ? + l). (8) 

For bipartite states of two qubits, p T A A B has at most one negative eigenvalue (Sanpera 
et al. 1998). For general multipartite and higher local dimension this is only a sufficient 
condition for the presence of entanglement. There exist entangled states with a positive 
partial transpose (PPT) known as bound entangled states (Acin et al. 2001, Horodecki 
et al. 1998a). The existence of bound entangled states ultimately limits the possibility 
to exploit the violation of Bell inequalities as a measure of entanglement. 

Let us briefly introduce the concept of bound entanglement. Such kind of 
entanglement can be recognized in terms of the so called distillation protocol, a non 
trivial procedure to optimize the extraction of Bell states from a mixture of entangled 
states (Bennett, Brassard, Popcscu, Schumacher, Smolin & Wootters 1996). The 
natural question, then is: Can any entangled state be actually distilled? The answer is 
yes for bipartite and qubit states (Horodecki et al. 1997). For multipartite entangled 
states and higher dimensional local Hilbert spaces a much more complex scenario 
emerges. In these cases examples of entangled states have been provided that cannot 
be distilled to maximally entangled states between the parties of the system, not even 
with an asymptotically infinite supply of copies of the state. Such a demoted form of 
entanglement was termed as bound entanglement. PPT entangled states were found 
first in (Horodecki et al. 1998 a). The existence of bipartite (with higher dimensional 
local Hilbert spaces) NPT bound entanglement has not excluded yet. This question 
has important implication on the additivity property of the distillable entanglement: 
if affirmative, entangled states could be generated from a mixture of non-distillablc 
states. Besides its speculative interest, it was demonstrated that bound entanglement 
can be activated in several quantum information and teleportation tasks to 'restore' 
the singlet fidelity of a given state (see (Horodecki et al. 2007)). 

Multipartite bound entangled states exist that are not fully separable but 
contain entanglement between each of their parties that cannot be distilled. Given 
that violation of PPT is necessary for distillation, a feature of multipartite bound 
entanglement is related to the 'incomplete separability' of the state (see (Dur & 
Cirac 2000)). A tripartite system A-B-C, for example, is separable with respect 
to the partition A — BC and B — AC and non-separable with respect to C — AB. 
The 'incomplete separability' is a sufficient condition for a state to have bound 
entanglement since the three qubits arc entangled and no maximally entangled 
state can be created between any of the parties by LOCC. For example, no 
entanglement can be distilled between C and A because no entanglement can be 
created with respect to the partition A — BC by LOCC. In the nest section the 
feature of incomplete separability will be exploited to detect multipartite bound 
entanglement in spin systems. The nature of correlations in bound entangled states 
is peculiar having both quantum and classical features. Therefore, if it is true that 
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all entangled states violates Bell inequality, the vice versa has not proven (see also 
(Popescu 1995) and Gisin (Gisin 1996)). Bipartite bound entangled states seem 
not violating Bell inequalities (Masanes 2006) (but some examples of multipartite 
bound entangled state violating Bell-type inequalities were found (Diir 2001)). For 
both the bipartite and multipartite entangled states there is a strong believe that 
necessary and sufficient condition for local realism is that the state satisfies the Peres 
criterium (Peres 1999, Werner & Wolf 2001, Werner & Wolf 2000). 

As already mentioned, a classification of multipartite entanglement is still 
missing. Nevertheless there are several quantities serving as indicators for multipartite 
entanglement when the whole system is in a pure state. The entropy of entanglement 
is an example for such a quantity and several works use multipartite measures 
constructed from and related to it (see e.g. (Coffman ct al. 2000, Meyer & Wallach 
2002, Barnum ct al. 2003, Scott 2004, de Oliveira, Rigolin & de Oliveira 2006, Love 
et al. 2007)). These measures give indication on a global correlation without discerning 
among the different entanglement classes encoded in the state of the system. The 
geometric measure of entanglement quantifies the entanglement of a pure state 
through the minimal distance of the state from the set of pure product states (Vedral 
et al. 1997, Wei & Goldbart 2003) 

E g {^) = -log 2 max | (*|$) | 2 (9) 

where the maximum is on all product states It is zero for separable states and rises 
up to unity for e.g. the maximally entangled n-particle GHZ states. The difficult task 
in its evaluation is the maximization over all possible separable states and of course the 
convex roof extension to mixed states. Despite these complications, a clever use of the 
symmetries of the problem renders this task accessible by substantially reducing the 
number of parameters over which the maximization has to be performed. Another 
example for the collective measures of multipartite are the measures introduced 
by Meyer and Wallach (Meyer & Wallach 2002) and by Barnum et al (Barnum 
et al. 2003, Barnum et al. 2004). In the case of qubit system the measure of Meyer and 
Wallach is the average purity (which is the average one-tangle in (Coffman ct al. 2000)) 
of the state (Meyer & Wallach 2002, Brcnnen 2003, Barnum ct al. 2004) 

2 N 

^ = 2 -^E T ^ 2 - ( 10 ) 

The notion of generalized entanglement introduced in (Barnum et al. 2003, Barnum 
et al. 2004) relaxes the typically chosen partition into local subsystems in real space. 
For the state \ip) it is defined as 

PA = Tr{[p A mm 2 } (ii) 

where Va is the projection map p — > Va(p). If the set of observables is defined by the 
operator basis {Ai,A 2 , . . . ,Al} then Pj, = J2f=i(Ai) 2 from which the reduction to 
Eq.(10) in the case of all local observables is evident. This conceptually corresponds 
to a redefinition of locality as induced by the distinguished observable set. Finally 
we mention the approach pursued in (Giihnc ct al. 2005) where different bounds 
on the average energy of a given system were obtained for different types of n- 
particlc quantum correlated states. A violation of these bounds then implies the 
presence of multipartite entanglement in the system. The starting point of Giihnc 
et al. is the notion of n- separability and k-producibility which admit to discriminate 
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particular types of n-particle correlations present in the system. A pure state | V) 
of a quantum systems of N parties is said to be n-separablc if it is possible to 
find a partition of the system for which | ip) = \(fii)\(p2) ■ ■ ■ \<t>n)- A pure state 
| ifi) can be produced by k-party entanglement (i.e. it is k-producible) if we can 
write | ip) = 101)1^2) ' ' ' \4>m) where the \<pi) are states of maximally k parties; by 
definition m > N/k. It implies that it is sufficient to generate specific k-party 
entanglement to construct the desired state. Both these indicators for multipartite 
entanglement are collective, since they are based on the property of a given many 
particle state to be factorized into smaller parts, k-separability and -producibility 
both do not distinguish between different k-particle entanglement classes (as e.g. the 
k-particlc W-states and different k-particlc graph states (Hcin ct al. 2004), like the 
GHZ state). Another approach is based on the observed relation between entanglement 
measures and SL(2,€>) invariant antilinear operators (Uhlmann 2000, Ostcrloh & 
Siewert 2005, Osterloh & Siewert 2006, Ostcrloh & Djokovic 2009). This allows 
certain sensitivity to different classes of multipartite entanglement (see also (Verstraete 
et al. 2002, Lamata et al. 2007, Bastin et al. 2009, Osterloh n.d.)). 

We close this section by reviewing how to swap or transmute different types of 
multipartite entanglement in a many body system into pairwise entanglement between 
two parties by means of generalized measures on the rest of the system. In a system of 
interacting spins on a lattice one could then try to maximize the entanglement between 
two spins (at positions i and j) by performing measurements on all the others. The 
system is then partitioned in three regions: the sites i, j and the rest of the lattice. This 
concentrated pairwise entanglement can then be used e.g. for quantum information 
processing. A standard example is that of a GHZ state (l/\/2)(| 000) + | 111)). After a 
projective measure in x-direction on one of the sites such a state is transformed into a 
Bell state. The concept of localizable entanglement has been introduced in (Verstraete, 
Martin-Delgado & Cirac 2004, Popp et al. 2005). It is defined as the maximal amount 
of entanglement that can be localized, on average, by doing local measurements in the 
rest of the system. In the case of N parties, the possible outcomes of the measurements 
on the remaining N — 2 particles are pure states \ip a ) with corresponding probabilities 
p s . The localizable entanglement E\ oc on the sites i and j is defined as the maximum 
of the average entanglement over all possible outcome states | ip s )ij 



where £ is the set of all possible outcomes (p s , \ips)) of the measurements, and E 
represents the chosen measure of entanglement of a pure state of two qubits (e.g. 
the concurrence). Although very difficult to compute, lower and upper bounds have 
been found which allow to deduce a number of non trivial properties of the state. 
An upper bound to the localizable entanglement is given by the entanglement of 
assistance (Laustsen et al. 2003) obtained from localizable entanglement when also 
global and joint measurements were allowed on the N — 2 spins . A lower bound of 
the localizable entanglement (Verstraete, Martin-Delgado & Cirac 2004) is fixed by 
the maximal correlation function between the two parties. 

4. Entanglement and magnetic order 




(12) 
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The entanglement present in the equilibrium (thermal or ground) state of a quantum 
system is very sensitive to the underlying collective behavior. This suggests, in the 
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case of spin systems, to analyze the relation between entanglement and magnetic 
order. We will discuss various aspects of this connection starting from the pairwise 
entanglement, we then proceed with the properties of multipartite entanglement. Most 
of the investigations available in the literature are for one-dimensional systems where 
exact results are available, later on we will overview the status in the d-dimensional 
case. 

The body of knowledge acquired so far makes it evident that entanglement in 
the ground state contains relevant information on zero temperature phase diagram of 
the system. We will highlight this relation in two paradigmatic cases when the spin 
system is either close to a quantum phase transitions or to factorizing field. We will 
mostly be concerned with a XYZ spin models in an external field. 

Range of pairwise entanglement As we discussed in Sec. 2, in an interacting spin 
system there exists a particular choice of the coupling constants and the external 
field for which the ground state is factorized (Kurmann et al. 1982, Giampaolo 
et al. 2009), i.e. the entanglement vanishes exactly. Several works were devoted 
to the characterization of the entanglement close to the factorizing point. It was 
demonstrated for one dimensional spin models in the class XXZ that the point 
at which the state of the system becomes separable marks an exchange of parallel 
and anti-parallel sector in the ground state concurrence (Fubini et al. 2006, Amico 
et al. 2006a). This change occurs through a global (long-range) reorganization of the 
state of the system. The range R of the concurrence (defined through the maximum 
distance between two sites over which it is non-zero) diverges. For the XY model it 
was found that this range is 



The existence of such a divergence has been confirmed in other one dimensional 
systems both for short (Amico et al. 2006a, Roscilde et al. 2004, Roscilde et al. 2005a) 
and long range interactions (Dusucl & Vidal 2005). This divergence suggests, as a 
consequence of the monogamy of the entanglement (Coffman et al. 2000, Osborne & 
Verstraete 2006), that the role of pairwise entanglement is enhanced on approaching 
the factorizing field (Roscilde et al. 2004, Roscilde et al. 2005a, Roscilde et al. 20056). 
Indeed, for the Ising model (i.e. 7=1), one finds that in this region the ratio between 
the two-tangle and the one-tangle tends to one(Fig.3) (Amico et al. 2006a). 

The diverging entanglement length is particularly intriguing in systems 
characterized by topological order. In Ref.(Son et al. 2009) the entanglement in the 
quasi-long-range ordered ground state of the one dimensional isotropic XY model. 
Because of the presence of the characteristic edge states in systems with non trivial 
topology, it was found that the quasi-long-range order is traced by entangled states 
localized at the edges of the system. The range of the concurrence diverges also 
close to the saturation field for a spin system with inverse-square interaction of 
the Haldane-Shastry type(Haldane 1988, Shastry 1988), interpolating in a sense 
between nearest-neighbor and fully connected graph interactions. In particular, in 
Ref.(Giuliano et al. 2009) it was shown that, in the absence of external magnetic field 
the Haldane-Shastry spin system only displays nearest neighbor entanglement; while, 
by increasing the magnetic field, the bipartite entanglement between spins at greater 
distance increases, up to a situation where the spin system saturates, entering a fully 
polarized phase described by a completely separable ground state. Given the special 
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Figure 3. Onc-tanglc Entanglement ratio for the XY model: the dashed line 
indicates the critical field. Upper panel: Value of 75/71 as h — > hi, as a function 

of h^ Y = v/l — 7 2 , for different values of 7; the Inset shows the "entanglement 
phase-diagram" in the h — ■y plane: the shaded (white) area corresponds to 
antiparallel (parallel) pairwise entanglement and the full line is the separable 
circle; the crosses indicate the position of the data shown in the main upper 
panel, while the dotted line corresponds to the x-axis used in the lower panel. 
Lower panel: 75/71 versus h, for 7 = 0.5. [From(Amico et al. 2006&)] 



role of the inverse-square interaction in one dimensional fractional statistics(Comtet & 
et al. 1999), this study is also relevant for the understanding of the interplay between 
the statistics and the entanglement. 

The range of concurrence does not diverge only at the factorizing field. There 
are one-dimensional spin systems where the pairwise entanglement has qualitative 
different features as a function of the distance between the sites. An example is 
the long-distance entanglement observed in (Campos Venuti, Degli Esposti Boschi 
& Roncaglia 2006). Given a measure of entanglement E(pij), Campos Venuti et al 
showed that it is possible that E(pij) ^ when \i — j\ — > 00 in the ground state. 
Long-distance entanglement can be realized in various one-dimensional models as in 
the dimerized frustrated Heisenberg models or in the AKLT model. For these two 
models the entanglement is highly non-uniform and it is mainly concentrated in the 
end-to-end pair of the chain. 

Pairwise entanglement and quantum phase transitions A great number of papers 
have been devoted to the study of entanglement close to quantum phase transition 
(QPT). Close to the quantum critical point the system is characterized by a diverging 
correlation length £ which is responsible for the singular behavior of different physical 
observables. The behavior of correlation functions however is not necessarily related to 
the behavior of entanglement. It is worth to stress that the study of entanglement close 
to quantum critical points does not provide new understanding to the scaling theory of 
quantum phase transitions. Rather it may be useful in a deeper characterization of the 
ground state wave function of the many-body system undergoing a phase transition. 
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Figure 4. The change in the ground state wave function in the critical region 
is analyzed considering the derivative of the nearest neighbor concurrence as a 
function of the reduced coupling strength. The curves correspond to different 
lattice sizes. On increasing the system size, the minimum gets more pronounced. 
Also the position of the minimum changes and tends as (see the left side inset) 
towards the critical point where for an infinite system a logarithmic divergence 
is present. The right hand side inset shows the behavior of the concurrence itself 
for an infinite system. [From(Osterloh et al. 2002)] 



In this respect it is important to explore, for instance, how the entanglement depends 
on the order of the transition, or what is the role of the range of the interaction to 
establish the entanglement in the ground state. We start by considering exclusively 
the properties of pairwise entanglement. 

Pairwise entanglement close to quantum phase transitions was originally analyzed 
in (Osborne & Nielsen 2002, Osterloh et al. 2002) for the Ising model in one dimension. 
Below we summarize their results in this specific case. The concurrence tends to zero 
for A ^ 1 and A « 1, the ground state of the system is fully polarized along the 
£-axes (z-axes). Moreover the concurrence is zero unless the two sites are at most 
next-nearest neighbors, we therefore discuss only the nearest neighbor concurrence 
C(l). The concurrence itself is a smooth function of the coupling with a maximum 
close to the critical point (but not related to any property of the phase transition). 
The critical properties of the ground state are captured by the derivatives of the 
concurrence as a function of A. The results are shown in Fig. 4. In the thermodynamic 
limit d\C{\) diverges on approaching the critical value as 

0aC(1)~ ^ln|A-A c | . (14) 

For finite system the precursors of the critical behavior can be analyzed by means 
of finite size scaling of the derivative of the concurrence. Similar results have been 
obtained for the XY universality class (Osterloh et al. 2002). Remarkably, although 
the concurrence describes short-range properties, nevertheless scaling behavior typical 
of continuous phase transition emerges. 

Over the last years the properties of pairwise entanglement were intensively 
studied. It was evidenced how it depends on the order of transition and on 
the universality class of the system. The bulk of results obtained so far can 
be summarized in a 'Ehrenfest classification scheme for entanglement', ultimately 
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arising because of the formal relation between the correlation functions and the 
entanglement. A way to put this observation on a quantitative ground is provided 
by a generalized Hohenberg-Kohn theorem (Wu et al. 2006). Accordingly, the ground 
state energy can be considered as a unique function of the expectation values of 
certain observables. These, in turn, can be related to (the various derivatives of) 
a given entanglement measure (Wu et al. 2004, Campos Venuti, Degli Esposti Boschi, 
Roncaglia & Scaramucci 2006). It was indeed shown that, given an entanglement 
measure M related to reduced density operators of the system, first order phase 
transition are associated to the anomalies of M while second order phase transitions 
correspond to a singular behavior of the derivatives of M. Also quasi-long range 
order is captured by the behavior of the pairwise entanglement (Gu ct al. 2003, Son 
et al. 2009). Other singularities like those noticed in the concurrence for models 
with three-spin interactions (Yang 2005), are due to the non-analyticity intrinsic in 
the definition of the concurrence as a minimum of two analytic functions and the 
constant zero. This was then explicitly shown for the quantum Ising, XXZ, and 
LMG models (Wu et al. 2006). For the Ising model, for example, the divergence of 
the first derivative of the concurrence is determined by the non-analytical behavior 
of (S x S x )(Wu et al. 2004). A relevant caveat to this approach is constituted by the 
uniaxial-LMG model in a transverse field (with h y = and 7 = 0) that displays a 
first order QPT for h x = 0. The concurrence is continuous at the transition since it 
does not depend on the discontinuous elements of the reduced density matrix (Vidal 
et al. 2004). The relation between entanglement and criticality was also studied in the 
spin-1 XXZ with single ion anisotropy. It was established that the critical anomalies 
in the entropy experienced at the Haldanc-largc- D (if an axial anisotropy D Ylii^i) 2 
is added to the Hamiltonian in Eq.(2)) transition fans out from the singularity of the 
local order parameter ((S z ) 2 ) (Campos Venuti, Degli Esposti Boschi, Roncaglia & 
Scaramucci 2006). Spontaneous symmetry breaking can influence the entanglement 
in the ground state. Below the critical field, the concurrence is enhanced by the parity 
symmetry breaking (Osterloh et al. 2006). Recently it was demonstrated that such 
enhancement is particularly pronounced for multipartite entanglement close to the 
symmetry breaking. This result constitute a further indication that multipartite, and 
not bipartite, entanglement plays the main role to establish long-range correlations at 
the critical points (dc Oliveira et al. 2008). 

In higher dimensions nearly all the results were obtained by means of numerical 
simulations. The concurrence was computed for the two dimensional quantum XY 
and XXZ models (Syljuasen 2003a). The calculations were based on Quantum Monte 
Carlo simulations (Sandvik & Kurkij 1991, Syljuasen & Sandvik 2002). Although 
the concurrence for the 2d models results to be qualitatively very similar to the 
one-dimensional case, it is much smaller in magnitude. The monogamy limits 
the entanglement shared among the number of neighbor sites. The ground state 
entanglement in two dimensional XYZ model were analyzed in (Roscilde ct al. 20056) 
by means of quantum Monte Carlo simulations. The divergence of the derivative of 
the concurrence at the continuous phase transition, observed in d = 1, was confirmed; 
also in this case the range of the pairwise entanglement extends only to few lattice 
sites. By studying the one- and the two tangle of the system, it was proved that 
the QPT is characterized by a cusp-minimum in the entanglement ratio ti/t 2 . The 
cusp is ultimately due to the discontinuity of the derivative of n . The minimum in 
the ratio T1/T2 signals that the enhanced role of the multipartite entanglement in the 
mechanism driving the phase transition. Moreover by looking at the entanglement it 
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was found that the ground state can be factorized at certain value of the magnetic 
field. The existence of the factorizing field in d = 2 was proved rigorously for any 
2d — XYZ model in a bipartite lattice. Unexpectedly enough the relation implying 
the factorization is very similar to that one found in d = 1 . 

Pairwise entanglement at finite temperature At finite temperature excitations 
participate to entanglement that can become non-monotonous on increasing 
temperature or magnetic field (Arnesen et al. 2001, Gunlycke et al. 2001). Concurrence 
for thermal states, was calculated in several situations (Osborne & Nielsen 2002, 
Wang 2002 a, Wang & Zanardi 2002, Tribedi & Bose 2006, Asoudch & Karimipour 
2004, Canosa & Rossignoli 2006, Canosa & Rossignoli 2005, Rigolin 2004, Zhang & 
Zhu 2006, Wang & Wang 2006, Zhang & Li 2005). At finite temperatures but close to a 
quantum critical points, quantum fluctuations are essential to describe the properties 
of the systems (Sachdev 1999). For illustration let us consider a one-dimensional 
quantum XY in an external magnetic field. Although such system cannot exhibits 
any phase transitions at finite temperature, the very existence of the quantum critical 
point is reflected in the crossover behavior at T ^ 0. According to the standard 
nomenclature, the renormalized-classical regime evolves into the quantum disordered 
phase through the so called quantum critical region (Sachdev 1999). In the T — h 
plane a U-shaped phase diagram emerges, characterized by the crossover temperature 
customarily defined as T cross = |A _1 — A^ 1 !. For T <C T cross the thermal Dc Broglic 
length is much smaller than the average spacing of the excitations; therefore the 
correlation functions factorize in two contributions coming from quantum and thermal 
fluctuations separately. The quantum critical region is characterized by T » T cross . 
Here we are in the first regime and the correlation functions do not factorize. In 
this regime the interplay between quantum and thermal effects is the dominant 
phenomenon affecting the physical behavior of the system. Thermal entanglement 
close to the critical point of the quantum XY models was recently studied by some 
of us (Amico & Patane 2007). In analogy with the zero temperature case it was 
shown that the entanglement sensitivity to thermal and to quantum fluctuations 
obeys universal T ^ 0-scaling laws. The crossover to the quantum disordered and 
renormalized classical regimes in the entanglement has been analyzed through the 
study of derivatives of the concurrence d\C and drC. The thermal entanglement 
results to be very rigid when the quantum critical regime is accessed from the 
renormalized classical and quantum disordered regions of the phase diagram; such 
a 'stiffness' is reflected in a maximum in drC at T ^ Tcross • The maximum in the 
derivatives of the concurrence seems a general feature of the entanglement in the 
crossover regime (see for example (Stauber & Guinea 2004, Staubcr & Guinea 2006)). 
Due to the vanishing of the gap at the quantum critical point, in the region T ^> T cross 
an arbitrarily small temperature is immediately effective in the system (see Fig. 5). 
From the analysis of the quantum mutual information it emerges that the contribution 
of the classical correlations is negligible in the crossover, thus providing the indication 
that such a phenomenon is driven solely by the thermal entanglement. It is interesting 
to study how the existence of the factorizing field hf affects the thermal pairwise 
entanglement (vanishing at zero temperature) . It results that the two-tangle r 2 is still 
vanishing in a region of the h — T plane fanning out from hf, therefore, if present, the 
entanglement in the region must be shared between three or more parties. In contrast 
to the analysis of the ground state, at finite temperature one cannot characterize the 
two separate phases of parallel and antiparallel entanglement. 
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Figure 5. The effect of temperature on the anomalies originated from the critical 
divergence of the field-derivative of C(R) can be measured by dT[d a C(R)]. The 
density plot corresponds to 7 = 1 and R = 1. T = T* and T = Tm are drawn as 
dashed and thick lines respectively. Maxima below T* are found at Tm = f3T croBB 
with /3 ~ 0.290 ± 0.005 and they arc independent of 7 and R; the crossover 
behavior is enclosed in between the two flexes of dT[d a C(R)] at T c i T c 2\ such 
values are fixed to: T cl = (0.170 ± 0.005)T C ro SS and T c2 = (0.442 ± 0.005)T cross 
and found to be independent of 7 and R. For T smaller than T c i dT[d a C(R)] ~ 0. 
Scaling properties are inherited in dj'[d a C{R)\ from d a C(R)] [From (Amico & 
Patane 2007)] . 



Thermal entanglement witnesses In some case it is hard to quantify the entanglement 
in a many-body system. Moreover it seems in general difficult to relate clear 
observables to some of the entanglement measures. If one relaxes the requirement 
of quantifying the entanglement and asks only to know if a state is entangled or not 
then in some important cases there is a very appealing answer in terms of the so called 
entanglement witness. Interestingly enough it was shown that entanglement witnesses 
in spin systems can be related to thermodynamic quantities (Toth 2005, Brukner 
& Vedral n.d., Wu et al. 2005, Hide et al. 2007). For the isotropic XXX or XX 
Heisenberg model, if the inequality is fulfilled (with U the internal energy, M z the 
magnetization) 

\U + h'M'\ 1 

~N\J\ 4' (15) 

then the system is in an entangled state. Once the internal energy and magnetization 
are calculated then it is possible to verify in which range of the parameters of the 
system and the external temperature entanglement is present. Most important is the 
fact that these types of inequalities can be verified experimentally. It should be stressed 
that the analysis based on the entanglement witness could be applied to any model 
for which we can successfully obtain the partition function. This feature is the main 
advantage of using thermodynamic witnesses approach to detecting entanglement. 
This method for determining entanglement in solids within the models of Heisenberg 
interaction is useful in the cases where other methods fail due to incomplete knowledge 
of the system. This is the case when only the eigenvalues but not eigenstates of 
the Hamiltonian are known (which is the most usual case in solid state physics) 
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and thus no measure of entanglement can be computed. Furthermore, in the cases 
where we lack the complete description of the systems one can approach the problem 
experimentally and determine the value of the thermodynamical entanglement witness 
by performing appropriate measurements. It is important to emphasize that any other 
thermodynamical function of state could be a suitable witness, such as the magnetic 
susceptibility or heat capacity (Wiesniak et al. 2005) 



Localizable entanglement The study of localizablc entanglement in spin chains 
allows to find a tighter connection between the scales over which entanglement and 
correlations decay (Verstraete, Martin-Delgado & Cirac 2004, Popp et al. 2005, Popp 
et al. 2006). One expects that the procedure of entangling distant sites by a set of 
local measurements will be less effective as the distance between the two particles 
increases thus leading to a definition of entanglement length For a translational 
invariant system £e can be defined in analogy of the standard correlation lenght 

e -i ,. , Eioc{\i-j\) 

^ = -|4h>oo 10g ■ (16) 

By definition the entanglement length cannot be smaller than the correlation length, 
S.E > £> therefore at a second order phase transition the localizable entanglement 
length diverges. In addition there may also appear "transition points" associated 
solely to a divergence in £e- In order to avoid misinterpretations, it must be 
stressed that the localizable "classical" two-point correlations then diverge as well. 
For the Ising model in a transverse field it can be shown that (Verstraete, Popp & 

Cirac 2004) max a = x , y , z |Q«| < E ioc {i - j) < \ £ ± ^where sf = (l ± (S?S?)) 2 - 

±<S?)) 2 and Q% = (S*S?) - {S?)(S?) . In this case, the lower bound 
is determined by the two-point correlation function in the x-direction. In the 
disordered phase (A < 1) the ground state possesses a small degree of entanglement 
and consequently its entanglement length is finite. The situation changes at the 
other side of the critical point. Here, although the correlation length is finite, 
the entanglement length is infinite as asymptotically the correlation tends to a 
finite values. The divergence of £e indicates that the ground state is a globally 
entangled state, supporting the general idea that multipartite entanglement is most 
relevant at the critical point (Osborne & Nielsen 2002, Roscilde et al. 2004). The 
properties of localizable entanglement were further investigated for a spin-1/2 XXZ- 
chain in (Jin & Korcpin 2004, Popp et al. 2005) as a function of the anisotropy 
parameter A and of an externally applied magnetic field h. The authors used exact 
results for correlation functions relying on the integrability of the models to find 
the required bounds. The presence of the anisotropy further increases the lower 
bound of the localizablc entanglement. At the Berezinskii-Kosterlitz-Thouless critical 
point (A = 1) the lower bound of the nearest neighbor localizable entanglement 
shows a kink (Popp et al. 2005). As pointed out by the authors this might have 
implications in the general understanding of the Berezinskii-Kostcrlitz-Thoulcss phase 
transitions where the ground state energy and its derivatives are continuous as well as 
the concurrence. The localizable entanglement in two-dimensional XX Z model was 
discussed as well (Syljuasen 2003 b) by means of quantum Monte Carlo simulations. 
A lower bound has been determined by studying the maximum correlation function 
which for A > — 1 is Q x , the long-range (power law) decay of the correlation implying 
a long ranged localizable entanglement. 



CONTENTS 



20 



For half-integer spins, gapped non-degenerate ground states are characteristic for 
systems in a disordered phase (consider paramagnets for example) . A finite gap in the 
excitation spectrum of the system in the thermodynamic limit makes the correlations 
decaying exponentially. This is the Lieb-Schultz-Mattis theorem, establishing that, 
under general hypothesis, the ground state of a spin system is cither unique and 
gapless or gapped and degenerate (Lieb et al. 1961) (see (Hastings 2004) for recent 
results). It was a surprise, when Haldane discovered that systems of integer spins 
can violate this theorem (Haldane 1983 a, Haldane 19836). Accordingly the long 
range order can be replaced by the so called hidden order of topological nature. 
This kind of order is established in the system because certain solitonic type of 
excitation become gapless for integer spins (Mikeska 1995). This suggests to investigate 
whether the entanglement in the ground state might play some role in establishing 
the hidden order characteristic for the Haldane phases. An aspect that might be 
relevant to this aim was recently addressed by studying the localizable entanglement 
in AKLT models (Verstraete, Martin-Delgado & Cirac 2004). The ground state of 
this class of models is of the valence bond type. For this case it was demonstrated 
that a singlet state made of two spins- 1/2 located at the ends of the chain can be 
always realized. This implies that the localizable entanglement is long ranged despite 
the exponentially decaying correlation (Verstraete, Martin-Delgado & Cirac 2004). 
Furthermore the localizable entanglement can be related to the string order parameter. 
The valence-bond-solid phase order was further studied by looking at the hidden order 
in chains with more complicated topology. The von Neumann entropy was studied 
in spin-1 XX Z model with biquadratic interaction and single ion anisotropy in (Gu 
et al. 2006, Wang et al. 2005) and in (Campos Venuti, Degli Esposti Boschi, Roncaglia 
& Scaramucci 2006). Some of the features of the corresponding phase diagram arc 
captured. The Haldane transitions exhibited in the phase diagrams are marked by 
anomalies in the Von Neumann entropy; its maximum at the isotropic point is not 
related to any critical phenomenon (the system is gapped around such a point), but it 
is due to the equi-probability of the three spin-1 states occurring at that point (Campos 
Venuti, Degli Esposti Boschi, Roncaglia & Scaramucci 2006). Since the Berezinskii- 
Kosterlitz-Thouless transition separating the XY from the Haldane or large- D phases 
connects a gapless with a gapped regime, it was speculated that an anomaly in the 
entanglement should highlight such transition (Gu et al. 2006). 

Multipartite entanglement Although pairwise entanglement allow to capture 
important properties of the phase diagram it was evidenced that the spin systems 
are most generically in multipartite entangled state(Wang 20026, Stelmachovic & 
Buzek 2004, BruBet al. 2005). Just to make an example we point out that the first 
excited state above a ferromagnetic ground state is a W-state that is a well-known state 
in quantum information with a multipartite entanglement. Despite its importance, 
a quantitative description of multipartite entanglement constitutes a challenging 
problem in the current research. In many-body physics multipartite entanglement has 
been studied resorting to 'global' measures that most often cannot distinguish different 
types of multipartite entanglement each other (sec however (Giihne et al. 2005)). A 
first way to estimate multipartite entanglement in spin system is provided by the 
entanglement ratio t 2 /ti as the amount of two spin relative to global entanglement. 
It is interesting to compare the behavior of such quantititics for quantum critical 
and factorizing points of spin models. In fact it emerged that t 2 /ti is small close 
to quantum critical points. In contrast the entanglement ratio approach to one close 
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Figure 6. One-tangle t\ and the sum of squared concurrences T2 as a function of 
the applied magnetic field in z-direction for the XY Z model with exchange along 
y: J y = 0.25 (in unit of exchange along z). Inset: contributions to the concurrence 
between j-th neighbors; full symbols stand for C (j), open symbols for C 11 (j)- 
The dashed line marks the critical field h c . [From (Roscilde et al. 2004)] 



to the factorizing point. Close to quantum critical points the entanglement ratio was 
calculated numerically for Id — XYZ (Roscilde et al. 2004). (Fig. 6). The entanglement 
ratio was calculated close to factorizing points for the quantum Ising model. 

To address multipartite entanglement directly several routes have been suggested. 
Of interest is the analysis based on geometric entanglement. Most of the works till 
now concentrated on critical systems. It was proved that geometric entanglement 
obeys an area law that, most probably, coincides with the well established area 
law of the von Neumann entropy (Botero et al. 2007, Orus 2008, Orus, Dusuel & 
Vidal 2008, Orus 2008, Shi et al. 2009). Such a result provides a further evidence that 
the universal behaviour of the block entanglement close to a critical point traces back 
to multipartite entanglement. 

Although some of the proposed measures rely on n-point correlation functions 
there is no clear evidence on whether this is a general feature of multipartite 
entanglement. As important example of measure of multipartite entanglement relying 
on two-point function was proposed in (de Oliveira, Rigolin & de Oliveira 2006, de 
Oliveira, Rigolin, de Oliveira & Miranda 2006, Somma et al. 2004) as 



where Pj.j+i is the reduced density matrix associated to the sites j and j + I. 
Similarly one can consider also three-body reduced density matrices and construct 
the corresponding global entanglement measure. Although the precise form has not 
been established yet, the global entanglement E^\ generalization of Eq.(17) should 
be related to the set of reduced n-qubits density operators. According to Oliveira and 
coworkers, then the hierarchy of E^ might provide a comprehensive description of 
entanglement in many-body systems already for moderate values of n. 

Global entanglement is very sensitive to the existence of QPTs. As a paradigmatic 
example the authors analyzed the phase diagram in the anisotropy-magnetic field 
plane. By extending an earlier approach developed in (Wu et al. 2004), de Oliveira 
et al. also showed how the non-analytic behavior of E^ is related to that of the 




N-l 



N 



(17) 
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ground state energy. Note that from Eq.(17) it is possible to define an entanglement 
length proportional to the correlation length £. This differs considerably from that 
one defined by the localizable entanglement (see Eq.(16)); the latter is always bounded 
from below by the correlation length and can even be divergent where £ is finite. 

As discussed in (Facchi et al. 2006a, Facchi et al. 2006&, Costantini et al. 2006) 
the analysis of the average purity might not be sufficient and the analysis of the 
distribution of the purity for different partitions could give additional information. 
Rather than measuring multipartite entanglement in terms of a single number, one 
characterizes it by using a whole function. One studies the distribution function of 
the purity (or other measures of entanglement) over all bipartitions of the system. 
If the distribution is sufficiently regular, its average and variance will constitute 
characteristic features of the global entanglement: the average will determine the 
"amount" of global entanglement in the system, while the variance will measure 
how such entanglement is distributed. A smaller variance will correspond to a larger 
insensitivity to the choice of the bipartition and, therefore, will be characteristic for 
different types of multipartite entanglement. 

In (Patane et al. 2007) multipartite entanglement is studied with the aim to 
shed light on how entanglement is shared in a many-body system. For a quantum 
XY model the simplest multiparticle entanglement of a subsystem made of three 
arbitrary spins of the chain is considered; then bipartite entanglement between a 
spin and the other two with respect to all possible bipartitions. It is found that the 
block of two spins may be entangled with the external spin, despite the latter is not 
entangled directly with any of the two spins separately (see Fig. 8). Hence the range 
of such spin/block entanglement may extend further than the spin-spin entanglement 
range. It is plausible that increasing the size of the subsystem considered will increase 
the range of the multipartite entanglement. For instance, the range of spin/block 
entanglement will increase if we consider a larger block. Hence, a single spin can 
be entangled with more distant partners, if one allows to cluster them into a large 
enough block. It would be intriguing to study how spin/block entanglement and, in 
general, block/block entanglement between subsystems, scale increasing the size of 
blocks, especially exploring the connection with quantum criticality. We remark that 
such analysis would be different with respect to the well known block entropy setting, 
since in that case one is interested in the block/rest-of-the-system entanglement. 

Bound entanglement We would like to conclude this brief description of the relation 
between entanglement and magnetic order by analyzing in which cases interacting spin 
system are in a bound entangled state (see section (3). We shall see that such kind 
of peculiar entangled states are generically 'engineered' by a many-body system at 
equilibrium as it occurs naturally in certain region of the phase diagram for the 'last' 
entangled states before the complete separability is reached. In this sense the bound 
entanglement bridges between quantum and classical correlations. Being the bound 
entanglement a form of demoted entanglement, it appears when quantum correlations 
get weaker. Bound entangled states were found in both the ground and thermal 
states of anisotropic XY modcls(Patanc et al. 2007, Ferraro, Cavalcanti, Garcia-Saez 
& Acin A 2008, Cavalcanti et al. 2008). We follow the approach pursued in (Patane 
et al. 2007) where three-spin entanglement in an infinite anisotropic XY chain was 
analyzed). Two different configurations were considered (see Fig. 7). 

At zero temperature bound entanglement appears (see Fig. 8) when the spins 
are sufficiently distant each others and as in the case of the spin/spin entanglement, 
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b) • • T 1 1 T T jtj t lT^ | T • • 

Figure 7. Configurations of spins described in the text and whose entanglement 
properties are presented in Fig. 8. We fix R = 3, hence a spin is directly entangled 
with its first three nearest neighbors, a) 'Clustering' two spins increases the range 
of entanglement (the scheme is symmetric also for spins on the left of the marked 
one), b) Symmetric configuration of spins such that no two- particle entanglement 
is present, but still the spins share multiparticle entanglement. 

it can be arbitrary long ranged near the factorizing field. To prove it, the idea is to 
resort the "incomplete separability" condition described in the first section. In fact 
from Fig. 8 we see that Msxt may be zero even if Ncentr is non-zero. Thus in such 
case the density matrix of the spins is PPT for the two symmetric bipartitions of one 
external spin vs the other two ("f | jf and jf | j) and Negative Partial Transpose 
(NPT) for the partition of the central spin vs the other two (we remark that PPT 
does not ensure the separability of the two partitions for dimensions of local Hilbert 
space greater than two). In fact if we could be able to distill a maximally entangled 
state between two spins then one of two previous PPT partitions would be NPT and 
this cannot occur since PPT is invariant under LOCC (Horodecki et al. 1998 b, Vidal 
& Werner 2002). Quantum states must be 'mixed enough' to be bound entangled. 
In the ground state a source mixing is the trace over the other spins of the chain. 
However if the spins are near enough the reduced entanglement is free. It results that 
the effect of the thermal mixing can drive the T = free entanglement to bound 
entanglement (see Fig. 9). This behavior shown for the Ising model is also found 
for the entire class of quantum XY Hamiltonians with generic values of anisotropy 
(the temperature at which the different types of entanglement are decrease with 7. 
NPT thermal bound entanglement was also found by (Ferraro, Cavalcanti, Garcia- 
Saez & Acin A 2008, Cavalcanti et al. 2008) resorting to block entropies. The idea of 
Acin and coworkers was to calculate the block entropy for two different bipartitions: 
one in which the two subsystems are made of contiguous block of spins (called 
'half-half partition); the other groups all the spins labeled, say, by even indices in 
one subsystem and the remaining ones in the other (called 'even-odd' partition). 
Because of the area law, the entanglement in the even-odd partition is more robust 
to thermal fluctuations than entanglement in the half-half partition (an hypothesis 
corroborated by actual calculations by the same authors). Therefore there is a range of 
temperatures for which the PPT condition is reached with an even-odd entanglement. 
The entanglement is bound because single particles cannot distill entanglement (as the 
half- half bipartion can be singled out to have particles in two different blocks). The 
calculations are done for a finite set of spins interacting according to an isotropic XY 
Hamiltonian. We remark that such NPT bound entanglement was found in closed 
systems. Recently it was demonstrated that NPT bound entanglement can arise 
dynamically from decoherence of multipartite entangled states of GHZ type(Aolita 
et al. 2008). Based on that it is intriguing to conjecture that such kind of entanglement 
could be generated dynamically via decoherence in open systems. Besides NPT, also 
PPT bound entangled states were found in spin systems at finite temperature(Toth 
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Figure 8. T = Negativities between one spin and the other two Vs magnetic 
field h are shown for both configurations of Fig 7. We consider 7 = 0.5. In this 
case outside the interval marked by the two solid vertical lines the range of spin- 
spin entanglement is R < 3 (for values of h inside this interval R grows due to the 
its divergence at factorizing field hf = yj\ — -y 2 ~ 0.86 (Amico et al. 2006a)). For 
configuration a) (upper panel), Nsiock signals the spin/block entanglement for a 
distance d = 4. For values of h outside the vertical lines, the Negativity signals 
genuine spin/block entanglement. For configuration b) (lower panel), both N~Ext 
(solid line) and Ncentr (dashed line) are plotted. For values of h outside the 
vertical lines the three spins share no spin/spin entanglement, hence for non zero 
N~Ext an d Ncentr f ree multiparticle entanglement is present. The latter turns in 
to bound entanglement for values of h such that only Afcentr 7^ and ATsxt = 
(both on the left and on the right of the solid lines). [From(Patane et al. 2007)] 



et al. 2007, Toth et al. 2009). The method developed by Toth et. al relies on certain 
relation between entanglement and squeezing of collective spins that can serve as 
separability test for separability of the given state (like an entanglement witness). 
The entanglement detected can be of multipartite type despite the relations involve 
only two-point correlation functions. They considered spin models of the Heisenberg 
type at finite size and proved that a range of temperature exists where their thermal 
state display multipartite entanglement that cannot be distilled for any bipartition of 
the system. 
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Figure 9. Threshold temperatures above that the negativity is zero in the even- 
odd (solid line), 1 : n — 1 (dotted line) and half-half (dashed line) partitions. 
We plot the threshold temperature as a function of the coupling parameter J of 
the model (1) with J^^ = Jy^^ = an d J^' 1 = OVij for n = 10 spins 

and /i z = 1.9. Inset: Temperature above which the negativity in the even- 
odd (solid line), 1 : n — 1 (dotted line) and half- half (dashed line) partitions is 
zero as a function of the number n of the spin Hamiltonian (1) with J = 1 and 
h z = 1.9.[From(Cavalcanti et al. 2008)] 



5. Conclusions and outlook 

The use of concepts developed in quantum information science has provided a new 
twist to the study of many-body systems. Here we presented a specific example of 
this kind of approach by discussing the relation between magnetism and entanglement. 
Looking at the next future it seems to us that the most challenging problems are a 
wider characterization of the multipartite entanglement and, in our opinion most 
important, a connection between this acquired knowledge and new experiments. 
Remarkable impact of quantum information in condensed matter has been proving 
on the possibility to design more efficient classical numerical algorithms for quantum 
many-body systems. 

As for experimental tests on entanglement in many-body systems, we observe that 
the most direct method seems relying on the entanglement witnesses that have been 
derived using thermodynamical quantities. Nevertheless methods based on neutron 
scattering techniques on magnetic compounds, that are of particular relevance for 
spin systems, are also valuable especially for a direct quantification of entanglement in 
macroscopic systems. In this context we notice that more refined experimental analysis 
seem to be required to extract entanglement. These might disclose new unexplored 
features of entangled many-body states. 
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